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Eaiianon (METHODS: In this article, the parameters of the Weibull distribution and the net 
Net premiuin premium are estimated using moment estimation, maximum likelihood, least squares 
Outliers of error, weighted least squares, percentage, Cramer-Von- Mises, Mixture of moment 
Weibull distribution 


and maximum-likelinood are discussed in the presence of outliers. R software was 
used for simulation and numerical calculations and Easyfit software was used to fit 
the Weibull distribution to the real example data. At the end, real data examples for 
obtaining different estimators of the premium in the case of unknown parameters B 
and 8 and known ais provided. 

FINDINGS: In this research bias, the mean square error of net premium and unknown 
parameters B and 9 were obtained by using different estimators for Weibull distribution 
data as well as the generalized variance of unknown parameters B and 8. 
CONCLUSION: In this part, the evaluation and comparison of the estimators using 
real and simulated data was discussed, which was obtained by different k_maxfor real 
data. For example, in the moment method, kmax was equal to 5, based on which the 
net premium is 3.37657. In the simulated data, according to k (number of outliers), 
n (sample size) and B and 9 values, bias values, mean squared error and generalized 
variance of premium and different estimators were obtained. As an example, for n=10, 
=1, B=1.5, 6=3 and a=70, by comparing the bias and generalized variance of the 
estimators, we come to the conclusion that based on the bias, the percentile estimator 
has a better performance than the other estimators. That is, it has less bias and 
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k=1@=70,8=15,0=3 lle BQ Mlered oilily \ Jone 
Table 1: Generalized variance of Bo and @ k=1, a=70, B=1.5, 8=3 


80 70 60 50 40 30 20 10 
3.87100E-05 1.26922E-06 6.06231E-07 3.21285E-07 2.14166E-07 1.69703E-07 1.55588E-07 1.92214E-07 GVinm 
1.97538E-08 1.99333E-08 2.46108E-08 2.84662E-08 3.65628E-08 5.22965E-08  8.36677E-08 — 1.79596E-07 GVnt 
3.11076E-08 3.15600E-08 3.85760E-08 4.25733E-08 5.63784E-08 7.93896E-08 1.21521E-07 2.44235E-07 GVinix 
7.97425E-09 5.61779E-06 4.42858E-05 7.15781E-06 2.35275E-05 1.68057E-03 4.80476E-02 1.23522E-01 GV; 
2.45829E-08 2.46707E-08 2.56103E-06 9.78616E-05 9.19961E-03 4.59410E-02 1.32637E-01 8.15719E-02 GVyis 
2.93741E-08 1.46580E-09 6.78105E-07 2.93950E-05 3.05403E-05 6.63587E-02 4.79328E-05 4.39627E-02 GV, 
1.91250E-09 5.42587E-08 2.81324E-08 3.09258E-08 1.11028E-08 1.26374E-05 5.50953E-05 3.15833E-04 GVeum 
8.01101E-06 6.95291E-06 8.25041E-06 8.91714E-06 9.42261E-06 9.42261E-06 9.24930E-06 9.55168E-06 GVaa 
k=1 5 @=70,B=15,0=3 lj O 5 B donde ashY Spur 
Table 2: Bias of premium, and @ for a= 70,8 =1.5,0 =3andk =1 
80 70 60 50 40 30 20 10 Shy 
6.613597 4.608409 3.402212 2.593288 2.013110 1.576279 1.235106 0.955721 Brn 
-0.020410 -0.021725 -0.023361 -0.025361 -0.028951 -0.034674 -0.045903 -0.080321 Bn 
0.009641 0.009457 0.009198 0.008803 0.008276 0.007367 0.005515 -0.000010 Bim 
0.013005 0.012189 0.012264 0.012404 0.012592 0.012583 0.012992 0.011409 Brat 
0.000300 0.000245 0.000301 0.000487 0.000497 0.000675 0.001127 0.002528 Oni 
0.009589 0.009397 0.009124 0.008731 0.008170 0.007218 0.005302 -0.000420 Brn 
0.023521 0.023343 0.024286 0.025545 0.027582 0.030612 0.037074 0.054289 Bes 
-0.020410 -0.021725 -0.023361 -0.025361 -0.028951 -0.034674 -0.045903 -0.080321 mix 
0.011859 0.011802 0.011709 0.011547 0.011365 0.011029 0.010302 0.008020 Prix 
0.000038 0.000299 0.000659 0.015561 0.029723 0.039065 0.056561 0.152185 Bis 
-0.003550 -0.003666 -0.007294 -0.002184 -0.006010 -0.003766 -0.013328 -0.014508 bis 
0.009245 0.009061 0.008408 0.008489 0.007509 0.006806 0.003836 -0.002298 Bis 
0.038261 0.039054 0.044906 0.061221 0.072476 0.106568 0.179438 0.164468 Buns 
0.000085 -0,000022 -0.000053 0.000231 -0.001159 -0.003675 -0.007245 -0.013306 Buts 
0.009568 0.009372 0.009097 0.008679 0.007986 0.006752 0.004361 -0.002463 Bwis 
0.000078 0.000094 0.000322 0.024120 0.042782 0.077545 0.057565 0.162843 B, 
-0.003718 -0,000092 -0.005916 0.000690 -0.002708 -0.000921 -0.004132 0.002380 6, 
0.009000 0.009195 0.008256 0.008415 0.007366 0.006419 0.003743 -0.002876 By 
0.000005 0.000015 0.000011 -0.000012 -0,004159 0.000115 0.003230 0.002263 Boom 
-0.002246 -0,004731 -0,008341 -0.004387 -0.002177 -0.005818 -0.012754 -0.011819 Bev 
0.009365 0.008914 0.008259 0.008249 0.007997 0.006678 0.004017 -0.001599 Bevin 
0.520921 0.508096 0.357149 0.430097 0.424165 0.428637 0.368719 0.403137 Baa 
-1.578636 -1,572513 -1,531119 “1.546045 -1.546783 -1,544725 -1.531580 -1,540284 Baa 
0.377298 0.373518 0.352403 0.358118 0.356877 0.353100 0.342251 0.329843 Pad 
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Table 3: Mean square error of insurance premium, B and 6 for a=70, B=1.5, 0=3 andk =1 


80 70 60 50 40 30 20 10 is ts al 
9' 
43.852004 21.270536 11.588148 6.730986 4.055707 2.486428 1.526584 0.914074 pee. 
0.000454 0.000513 0.000594 0.000700 0.000909 0.001299 0.002249 0.006739 On 
0.000093 0.000090 0.000085 0.000078 0.000069 0.000055 0.000032 0.000003 Bmm 
0.000988 0.000928 0.000943 0.000900 0.000952 0.000980 0.001033 0.001006 Brat 
0.000024 0.000026 0.000032 0.000039 0.000047 0.000066 0.000101 0.000226 Bat 
0.000092 0.000089 0.000084 0.000077 0.000067 0.000053 0.000029 0.000002 Pmt 
0.001378 0.001330 0.001389 0.001402 0.001557 0.001763 0.002239 0.003824 Bax 
0.000454 0.000513 0.000594 0.000700 0.000909 0.001299 0.002249 0.006739 Bints 
0.000141 0.000140 0.000138 0.000134 0.000130 0.000123 0.000108 0.000068 Dmix 
0.000003 0.000408 0.001299 0.000739 0.001538 0.079502 0.612501 1.460405 Bis 
0.021561 0.021594 0.043122 0.014423 0.035980 0.021606 0.079039 0.086155 61, 
0.000347 0.000342 0.000590 0.000246 0.000488 0.000304 0.000945 0.000994 Bis 
0.002323 0.002343 0.070741 2.139438 1.274108 2.135238 3.111170 1.061570 Burts 
0.000029 0.000030 0.000037 0.000046 0.007252 0.021644 0.043131 0.079042 bvis 
0.000092 0.000088 0.000083 0.000076 0.000150 0.000303 0.000529 0.000971 Bwis 
0.000007 0.000043 0.000024 0.004645 0.002897 3.074393 0.004270 1.039910 By 
0.028742 0.007230 0.043148 0.007235 0.028801 0.021628 0.050318 0.043388 6, 
0.000429 0.000171 0.000588 0.000158 0.000399 0.000299 0.000607 0.000503 Bp 
0.000000 0.000002 0.000001 0.000001 0.000094 0.000352 0.000710 0.004011 Boom 
0.014401 0.028755 0.050319 0.028763 0.014417 0.035954 0.079013 0.078996 Osim 
0.000262 0.000427 0.000675 0.000415 0.000237 0.000475 0.000949 0.000909 Devin 
0.300706 0.286001 0.139618 0.200382 0.194808 0.199080 0.149837 0.177763 Baa 
2.501719 2.481907 2.350353 2.396420 2.398889 2.392302 2.352259 2.379324 Ona 
0.144818 0.141826 0.125477 0.129700 0.128803 0.126068 0.118449 0.110076 Bad 
k=3 5 @=70,B=15,0=3 hla O 5 B dlimrasd willy F Jr 
Table 4: Generalized variance of pand 6 for @=70, 8 =1.5,0 =3and k =3 
n 
80 70 60 50 40 30 20 10 
lasj3l 

-1.0912E08 -1.897E08  -1.9171E08 -1.5905E08 1.6600E-08  2.0900EF-08  2.6200E-08  4.4200E-08 GVinm 

-5.2187E09 -7.4700E09  -7.5899E09 = -7.3745E09 = 8.4200E-09 1.3900E-08  1.9300E-08  5.3100E-08 GVini 

-1.0017E08 -1.0187E08  -1.1286E08  -1.0466E08 1.2600E-08 1.9500E-08  2.7200E-08 5.5000E-08 GV nix 

-0.0817E08 -1.1317E08  -1.4030E08 -2.6497E05 5.1500E-06  2.4500E-08 4.9900E-05  5.4200E-05 GVis 

-1.0313E08 -1.1613E08  -1.2515E08  -1.3969E08 1.7200E-08  8.6000E-05 5.4400E-06  5.6300E-06 GVyis 

-1.1153E08 -1.1456E08  =-1.1821E08  -1.0064E05 5.7900E-03 5.9181E-03 1.0552E-04  8.6500E-06 GY, 

-2.1246E09 -1.0146E08  -4.9640E05 -4.6588E06 1.3800E-05 7.5900E-06 9.7800E-06  9.3300E-08 GVveum 

-1.0013E05 -1.0878E05 -1.0034E05_—s--7.1534E05 ~—7.4300E-06  8.5700E-06  9.1200E-06 7.3300E-06 GVad 
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Table 5: Bias of premium B and 6 for g 70, B =1.5,0 =3and k =3 


n 
80 70 60 50 40 30 20 10 
yaly 
1.470136 1.342069 1.234950 1.136308 1.043646 0.955707 0.870904 0.778491 Bm 
-0.039010 -0.04101 -0.04544 -0.053053 -0.063128 -0.080208 -0.115863 -0.227768 Bnm 
0.027328 0.028338 0.027505 0.026485 0.024802 -0.000030 -0.005519 -0.022033 Dini 
0.0056801 0.004580 0.003757 0.003776 0.003150 0.002876 0.002904 0.001890 Brn 
0.000132 0.000503 0.000985 0.000281 0.000525 0.001172 0.001311 0.003346 Onl 
0.029350 0.028260 0.027422 0.026391 0.024711 -0.000163 -0.005706 -0.022443 Bm 
0.023550 0.025690 0.027387 0.030992 0.035729 0.044763 0.063962 0.127084 Brix 
-0.039012 -0.04101 -0.04544 -0.053053 -0.063128 -0.080208 -0.115863 -0.227768 O nix 
0.033010 0.032700 0.032358 0.032016 0.031339 0.008130 0.005710 -0.003167 Pmix 
0.015698 0.014598 0.013217 0.012882 0.012238 0.011512 0.011847 0.011612 Bis 
-0.000105 -0.00020 0.000413 -0.015268 -0.003157 0.000107 -0.015360 -0.011244 an 
0.031235 0.028275 0.027416 0.024753 0.024364 -0.000174 -0.007286 -0.023127 Dis 
0.015267 0.014258 0.012579 0.012096 0.011509 0.009486 0.010852 0.010747 Beis 
0.0000520 0.000022 0.000534 -0.000191 -0.000018 -0.002696 -0.003207 -0.004972 wis 
0.029147 0.028253 0.027408 0.026370 0.024679 -0.000525 -0.006025 -0.022543 Bwis 
0.015311 0.014612 0.013145 0.012780 0.028513 0.033665 0.014419 0.009138 By 
0.0009523 0.001034 0.001872 -0.004512 -0.000852 -0.002633 0.001883 0.009463 6, 
0.029061 0.028141 0.027259 0.025535 0.024100 -0.001222 -0.006604 -0.023883 Bp 
0.003852 0.002779 0.001828 0.002216 0.001621 0.000954 0.001558 0.001181 Boom 
-0.000052 -0.00004 -0.02927 -0.003062 -0.008789 -0.002352 -0.005404 0.001550 Boom 
0.032144 0.028331 0.024214 0.026119 0.023796 -0.000395 -0.006226 -0.022229 Bevm 
0.613400 0.512420 0.413913 0.344350 0.457913 0.440960 0.403016 0.349001 Baa 
-1.593504 -1.57290 -1.54372 -1.532222 1.549700 -1.546344 -1.542829 -1.535349 Gua 
0.386591 0.385821 0.369460 0.361285 0.363760 0.331751 0.314774 0.265815 Bad 
k=3 5 @=70,B=15,0=3 Gls O59 B dome GUS 999 lpi GeSilee ¥ Spar 
Table 6: Mean square error of insurance premium, B and © for g =70,8 =1.5,0 =3andk =3 
n 
80 70 60 50 40 30 20 10 
laaysi 
2.01253 1.801670 1.525489 1.291486 1.089438 0.913589 0.758654 0.606224 Bm 
0.001360 0.001730 0.002116 0.002869 0.004055 0.006531 0.013572 0.052149 Gam 
0.000856 0.000804 0.000757 0.000702 0.000616 0.000001 0.000032 0.000488 PBmm 
0.000275 0.000240 0.000244 0.000210 0.000197 0.000220 0.000210 0.000262 Brn 
0.000027 0.000031 0.000035 0.000040 0.000048 0.000075 0.000112 0.000306 Ont 
0.000851 0.000799 0.000752 0.000697 0.000611 0.000001 0.000034 0.000506 Pm 
0.000761 0.000881 0.000977 0.001158 0.001466 0.002216 0.004294 0.016399 Breve 
0.001456 0.001730 0.002116 0.002869 0.004055 0.006531 0.013572 0.052149 8B mnix 
0.001092 0.001070 0.001048 0.001026 0.000983 0.000067 0.000034 0.000013 Bmix 
0.000538 0.000528 0.000512 0.000463 0.000435 0.000424 0.000703 0.001978 Bis 
0.000032 0.000038 0.000042 0.089732 0.018042 0.000094 0.089638 0.053781 6, 
0.000852 0.000800 0.000752 0.001664 0.000803 0.000001 0.001042 0.001042 Dis 
0.000652 0.000558 0.000480 0.000462 0.000446 0.004886 0.000486 0.001877 Buis 
0.000029 0.000034 0.000039 0.000046 0.000058 0.018033 0.018027 0.018162 Bwis 
0.000852 0.000799 0.000752 0.000696 0.000610 0.000205 0.000232 0.000669 Bwis 
0.000588 0.000529 0.000456 0.000444 0.321345 0.166055 0.006054 0.002588 By 
0.000031 0.000039 0.000046 0.035922 0.018052 0.035891 0.018055 0.018468 6, 
0.000840 0.000792 0.000744 0.001074 0.000790 0.000414 0.000244 0.000729 By 
0.000301 0.000291 0.000281 0.000265 0.000260 0.000421 0.000276 0.000346 Bevm 
0.000052 0.000038 0.180193 0.017954 0.053900 0.018066 0.035930 0.000380 Boom 
0.000769 0.000803 0.002709 0.000894 0.001192 0.000206 0.000437 0.000497 Die 
0.351410 0.291318 0.185775 0.128892 0.231765 0.211974 0.177168 0.132000 Baa 
2.582651 2.483567 2.389717 2.352914 2.408907 2.397471 2.387199 2.362506 6aa 
0.164085 0.151088 0.137909 0.131525 0.133899 0.111312 0.100208 0.071215 Bad 
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Fig. 1: Bias, mean squared error and generalized variance of different estimators based on simulation results for 


B=1.5,0=3,a =70andk =1 
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Fig. 2: Bias, mean squared error and generalized variance of different estimators based on simulation results for 


B=1.5,0 =3,a =70and k =3 
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Fig. 3: Histogram, box plot, quantile-quantile and empirical distribution function of Danish fire damage data 
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Table 7: Various estimators of B for Q =1.3248 in respect to k 


Baa = By Burs Bis Bintx But Binm k 
2.793501 0.9999996 0.9999997 1.00002044 0.9999995 3.779443 =0.1336535 0.1254932 1 
2.793501 1.0000236 1.0000209 0.06839348 1.0000285 3.761079 2.8841202 0.2699722 2 
2.793501 1.0000243 1.0000214 1.00002152 1.0000293 3.742452 0.1195378 0.4379768 3 
2.793501 1.0000247 1.0000183 0.99999969 1.0000297 3.723491 2.9819105 0.6356001 4 
2.793501 1.0000255 1.0000225 1.00002267 1.0000305 3.704073 3.0357152 0.8712122 5 

Ke ay cnn 11.3248 ly Q cdlzee la Sojol yA Jose 
Table 8: Different estimators of 0 for @ =1.3248 respect to k 

Baa Bev 6, Owis 615 8 mix Ont 8 nm k 
6.32406 0.4112209 0.4274100 0.3614717 0.4120285 0.1910469 0.3219974 0.1910469 1 
6.32406 0.4112198 0.4274089 0.3873936 0.4120272 0.1870536 0.2344257 0.1870536 2 
6.32406 0.4112194 0.4274086 0.3614712 0.4120268 0.1830430 0.4966554 0.1830430 3 
6.32406 0.4112191 0.4274084 0.3614725 0.4120264 0.1790135 0.2247280 0.1790135 4 
6.32406 0.4112187 0.4274079 0.3614706 0.4120259 0.1749633 0.2197206 0.1749633 5 

k Cons ers cqilsicw yes ali 4 Jane 
Table 9: Empirical likelihood function with respect to k 
laa, Baa) U6 cam Bevm) 16, B») UGyis, Bus) 1(6,,, Bis) U6 mix) Bix) Unt Bri) Umm Bimm) k 
-2183.302 -76.22000 -77 83384 -71.92480 -76.29815 =67.79172 -61.97735 -67.69032 1 
-2215.391 -73.54580 -75.15964 -55.36909 -73.62395 -64.68653 -63.55078 -65.28727 2 
-2248.454 -71.31221 -72.92605 -67.01706 -71.39036 -62.01014 -50.28831 -63.75286 3 
-2282.949 -69.40267 -71.01651 -65.10756 -69.48082 -59.64774 -58.57895 -62.58078 4 
-2318.622 -67.75401 -69.36785 -63.45886 -67.83216 -57.53808 -56.50286 -61.31777 5 
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Table 10: Values of B... 8 esKonax and net premium of different estimators 
Kmax — Bmax = Boiss p 
kmm = 5 Bmm = 0.8712122 Imm = 0.1749633 Bmm = 3.37657 
km = 3 Bmi = 0.1195378 8 mi = 0.4966554 Bm = 3.217375 
kmix = 5 Brix = 3.704073 mix = 0.1749633 Binix = 4.619323 
kis =5 Bis = 1.000031 1 = 04120259 Bis = 2.233271 
Kwis = 2 Bwis = 0.06839348 Bwis = 0.3873936 Pwis = 4.532246 
ky =5 Pp = 1.000023 4 = 0.4274079 Bp = 2.152901 
Resa = Boom = 1.0000255 Deum = 0.4112187 Bevm = 2.237657 
kaa =1 Baa = 2.793501 aq = 6.32406 Baa = 0.1423892 
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Fig. 4: Histogram, box plot, quantile-quantile and empirical distribution function of Iran Insurance car insurance claim data 
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AS 95g) Codd; Alize clagig, 41) dan > Soyal ne oleic 
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S83) 3 Km Gly KH popeSle polio 1+ Jae 59 
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k=5 oolidljgs claools stax ogy G59) jo «Jie sgba 
k=S« 8 =0.8712122 chee Galply cul oad oygly 
Ol & bles Gd dag ge> 9 Cul 6 =0.1749633 5 
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Table 11: Various estimators 


of £ for a =1.372 respect to k 


Baa Bam By Bwis Bis Braise Bt cain k 
5.000085 4.9999666 4.9999742 4.999959 4.999990 3.142297 18.353170 0.4316495 1 
5.000085 0.1239608 1.5993283 1.780974 2.124528 3.099228 9.098362 1.1095494 2 
5.000085 4.9984942 0.2618328 1.819728 2.150974 3.055664 8.250227 2.3236309 3 
5.499902 4.9984942 1.6521680 5.000005 2.187310 3.008829 7.195901 5.1099639 4 
5.499902 5.0000020 1.6891253 4.999372 5.000003 2.960044 7.813507 17.9738177 5 

Ke ay Sacked OA NVY ely O Blede ele Faigle 2V¥ lace 
Table 12: Different estimators of 9 for @ =1.372 with respect to k 
Goa Bom 6, Owais brs 8 nix 8rnt Briss k 
4.355500E-07 4.676806E-07 —3.957359E-07 4.363649E-07 3.098155E-07 2.796704E-07  3.098155E-07. 1 
1.376386E-05 5.291053E-07  4.771573E-07 — 4.052705E-07 4.327867E-07 3.007932E-07 5.062374E-07 3.007932E-07 2 
1.376386E-05  3.680291E-07 5.825220E-07 —3.910359E-07 4.143691E-07 2.917207E-07 4.106268E-07 2.917207E-07. 3 
1.376386E-05  3.345015E-07 4.507850E-07 —2.941704E-07 3.958547E-07 2.825896E-07 4.777218E-07 2.825896E-07 4 
1.376386E-05 _3.026454E-07 _4.362148E-07 __2.655536E-07 3.025644E-07 2.733899E-07 1.541216E-07 07 5 
A aed ay alice ye QUAY Joue 
Table 13: Different estimators of 8 for a=1.372 with respect to k 
U(Gaas Baa) U(Bcom Bevm) 1(6,, Bo) L(8wis, Bwis) 1(615, Bis) UGrnixs Brix) L(G Bm) U(Grm> Bm) k 
-3369.084 -312.8188 -313.7861 -311.9326 -312.8407 -311.2676 -315.5855 -311.7098 1 
-3566.087 -307.9107 =311,0299 -309.2609 -309.8204 -308.7950 -318.8814 -309.4394 Z 
-3766.795 -308.2520 -305.2523 -307.1960 -307.6209 -306.7795 -313.5469 -306.8922 3 
-Inf -306.4316 -307.1504 -305.8656 -305.7868 ~305,1112 -316.8777 -305.8934 4 
-Inf -304.8807 -305.6489 -304.5234 -304.8793 -303.7345 -308.8130 -324.9688 5 
bre Gl 505] » valle dan b> 2. Br? nas Kan polie VF Joe 
Table 14: Values of Bins Bo Knax and net premium of different estimators 
Kmax x Brvas — Brice p 
kimm = 4 Bm = 5.109964 nm = 2-825896E — 07 Pmm = 2462050 
km =5 Bt = 7.813507 6,1 = 1.541216E — 07 Pmt = 4639688 
Kmix = 5 Bmix = 2.960044 O mix = 2.733899E — 07 Bmix = 2791065 
kig=5 Bis = 5.000003 6,, = 3.025644E — 07 Bis = 2417798 
Rwis = 5 Buns = 4.999372 Oyis = 2.655536E — 07 Bwis = 2754793 
kp = 3 By = 0.2618328 6, = 5.82522E — 07 By = 2233599 
Kevm = 5 Boum = 5.000002 Bum = 3.026454E — 07 Bevm = 2417151 
Kaa = 1 Bag = 5.000085 64q = 1.376386E — 05 Baa = 63779.17 
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